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TEICHNICAL  NOTE  2903 

IMPINGEMENT  OP  CLOUD  DROPLETS  ON  AERODYNAMIC  BODIES  AS  AFFECTED 
BY  COMPRESSIBILITY  OF  AIR  FLOW  AROUND  THE  BODY 
By  Rinaldo  J.  Brun,  John  S.  Serafinl,  and  Helen  M.  GeLllagher 


SUMMARY 

The  trajectories  of  water  droplets  In  a compressihle-alr  flow  field 
around  a cylinder  were  computed  with  a mechanical  analog.  The  results  of 
the  calculations  at  approximately  the  flight  critical  Mach  number  were 
compared  with  calculations  of  trajectories  in  an  incon5>ressible  flow 
field.  For  a cylinder,  the  effect  of  conpressibility  of  the  air  on  the 
droplet  trajectories  was  negligible  up  to  the  flight  critical  Mach  number . 
The  results  obtained  with  the  cylinder  were  extended  to  airfoils.  This 
extension  is  possible  because  the  incompressible  flow  fields  of  both 
cylinders  and  airfoils  are  similarly  altered  by  compressibility. 


INTRODUCTION 

The  analyticsQ.  study  of  ice  formation  on  airfoils  and  other  aero- 
dynamic bodies  exposed  to  icing  clouds  requires  the  csuLctilation  of  the 
trajectories  of  the  droplets  with  respect  to  the  body.  Previous  inves- 
tigators have  calculated  the  water -droplet  trajectories  with  respect  to 
cylinders  (references  1 to  3)  and  with  respect  to  airfoils  (refer- 
ences 4 to  6).  All  the  calcvilations  in  the  references  cited  were  made 
with  the  assumption  that  the  air  surrounding  the  aerodynamic  bodies 
behaved  as  an  incompressible  fluid.  An  evaluation  of  the  effect  of 
the  compressibility  of  air  in  the  flow  field  on  the  trajectories  of 
cloud  droplets  is  important,  because  many  airfoils  are  now  operated  at 
high  speed. 

A knowledge  of  the  air  flow  around  the  aerodynamic  body  is  required 
in  the  calculations  of  the  droplet  trajectories,  because  the  motion  of 
a droplet  is  governed  by  the  drag  forces  imposed  on  the  droplet  by  its 
motion  relative  to  the  air  moving  along  the  air  streamlines.  The  droplet 
momentum  tends  to  keep  the  droplet  moving  in  a straight  path,  while  the 
drag  forces  tend  to  force  the  droplet  to  follow  the  streamlines. 
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Two-dinffinsioned  incon5>res8lble  flow  fields  may  be  determined  either 
from  known  mathematical  relations  or  from  electrolytic  tank  methods  for 
most  aerodynamic  bodies.  The  compressible-air  flow  fields  are  extremely 
difficult  to  obtain  for  most  aerodynamic  bodies  and  are  usually  aijproxi- 
matlons  or  incompressible  flow  fields  with  corrections  applied.  A 
reasonably  simple  method  for  obtaining  droplet- trajectory  results  in  a 
compressible-air  flow  field  is  to  evaluate  the  compressibility  correc- 
tion to  be  applied  to  trajectory  data  csQ-culated  for  an  inconqpresslble 
flow  field.  Although  the  compressible -air  flow  field  around  a cylinder, 
which  is  obtained  by  applying  correction  terms  to  the  incompressible 
flow  field,  is  laborious  to  calculate,  the  cylinder  flow  field  provides 
a relatively  simple  method  for  eveLLuatlng  the  correction  required  for 
compressibility. 

In  order  to  evaluate  the  magnitude  of  the  compressibility  effect 
on  droplet  impingement,  trajectories  in  a compressible  flow  field  around 
a cylinder  were  calculated  at  the  NACA  Lewis  laboratory  and  the  results 
compared  with  those  of  trajectories  calculated  in  an  incompressible  flow 
field  (reference  3).  The  compressible-air  flow  field  was  determined 
only  for  the  cylinder  flight  critical  Mach  number,  irtiich  is  defined  as 
that  airspeed  which  results  in  sonic  velocity  at  some  location  on  the 
cylinder . 

The  results  of  the  study  of  the  effect  of  the  compressibility  of 
air  around  a cylinder  are  extended  to  incl\ide  airfoils.  The  extension 
is  based  on  a comparison  of  the  compresslhle-air  flow  field  around  a 
cylinder  at  the  cylinder  flight  critical  Mach  numiber  with  the 
compressible-air  flow  field  around  an  airfoil  at  the  airfoil  critical 
Mach  numiber. 


SYMBOLS 

The  following  symibols  are  used  in  this  report: 

a droplet  radius,  ft  (3.048Xip^  microns)  

Cj3  drag  coefficient  for  droplets,  dimensionless 

E collection  efficiency  based  on  cylinder  radius,  dimensionless 

2 

K inertia  paraioeter,  ^ dimensionless 

a HqL 

L cylinder  radius  or  airfoil  chord  length,  ft 
M Mach  numiber,  dimensionless 


2450 


MCA  TS  2903 


3 


Re 


Reo 


r 


t 

U 

u 


V 


X. 


a 

r 

e 

n 

p 

T 


<? 


local  Reynolds  n\mber  with  respect  to  droplet,  2a  v/p,  dimension- 
less 


free-stream  Reynolds  number  with  respect  to  droplet, 
dimensionless 


radial  distance  from  cylinder  axis  to  point  in  flow  field,  ratio  to 
cylinder  radixis 

airfoil  thickness,  ratio  to  chord  length 
free-stream  velocity,  ft/ sec 

local  air  velocity,  ratio  of  the  actual  velocity  to  the  free-stream 
velocity 

local  droplet  velocity,  ratio  of  actual  velocity  to  free-stream 
velocity 

local  vector  difference  between  velocity  of  droplet  and  velocity 
of  fl\iid,  ft/ sec 

rectangular  coordinates,  ratio  of  actual  distance  to  cylinder 
radius  L or  airfoil  chord  length 

angle  between  x-axls  and  radial  line  to  point  in  flow  field,  deg 

ratio  of  specific  heats 

impingement  angle  on  cylinder,  deg 

viscosity  of  air,  slugs/(ft) (sec) 

density,  slugs/ cu  ft 

dimensionless  time_  scale,  ratio  to  time  required  for  droplet  to 
travel  distance  L at  velocity  U 


^^0^  = 18LU 
K - \ /o' 


dimensionless 


Subscripts: 
a air 


c coinpresslble 
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i incompressible 

m maximum  < 

s s\jrface 

w water 

X horizontal  component 
y vertical  con^ranent 
0 free  stream 


AMLYSIS 


The  differential  equations  that  describe  the  motion  of  a droplet 
with  respect  to  a cylinder  have  been  derived  in  reference  3 and  are  pre- 
sented here  in  the  following  form: 


dT 


Cl>Re  1 

“2r  K K - ^x> 


CpRe  1 
24  K 


(Uy  - 


(1) 


where 


K 3 


2 Pw^^^ 
5 jiqL 


(2) 


and  the  Reynolds  number  Re  can  be  obtained  conveniently  in  terms  of 
the  free -stream  Reynolds  number 


Rsq  = 2aU 


(3) 


such  that 


(4) 


The  differential  equations  (l)  state  that  the  motion  of  a droplet 
is  governed  by  the  drag  forces  imposed  on  the  droplet  by  virtue  of  the 
relative  velocities  between  the  droplet  and  the  local  air  in  motion  along 
the  streamlines  arotand  the  cylinder.  The  conpresslble-alr  flow  around  a 
right  circular  cylinder  of  infinite  extent  is  described  by  the  following 
equations : 
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The  first  bracketed  term  describes  the  incompressible  flow  around  the 

cylinder.  The  term  with  the  factor  and  the  term  with  the  M*  fac- 
tor are,  respectively,  the  first-  and  second-order  corrections  for  the 
effect  of  con^jressibility  on  the  Incoatpresslble  flow  field.  The  expres- 
sion for  the  velocity  potential,  from  which  equations  (5a)  and  (5b)  were 
derived,  was  obtained  from  reference  7.  In  equation  (5),  a.  is  measured 
counterclockwise  from  the  x-axis  of  the  cylinder,  as  is  shown  in  figure  1. 

Equations  (l)  to  (5)  are  written  in  dimensionless  form  in  order  to 
maintain  _the  number  of  calculations  to  a Tn^n^Tm1Tn  and  to  sin^llfy  the 
presentation  of  the  results.  The  equations  apply  to  the  motion  of  drop- 
lets in  a plane  perpendicular  to  the  axis  of  the  cylinder,  which  is 
located  at  the  origin  of  the  coordinate  system  (fig.  1).  At  an  infinite 
distance  ahead  of  the  cylinder  the  uniform  air  flow  carrying  the  cloud 
droplets  is  assumed  to  be  approaching  the  cylinder  from  the  negative 
x-direction  and  parallel  to  the  x-axis.  All  the  distances  appearing  in 
the  equations  are  ratios  to  the  cyliad.er  radius  L,  which  is  assumed  to 
be  the  unit  of  distance.  The  velocities  appear  as  fractional  parts  of 
the  free-stream  velocity  U.  Time  is  expressed  in  terms  of  the  cylinder 
radivis  and  free-stream  velocity  such  that  the  \anit  of  time  is  the  time 
required  for  a droplet  to  travel  a distance  L at  velocity  U. 

Tte  more  important _ assumptions  that  have^been  necessary  in  order  to 
solve  the  problems  are : 

(1)  At  a large  distance  ahead  of  the  airfoil  (free-stream  conditions) 
the  droplets  move  with  the  same  velocity  as  the  air. 

(2)  The  droplets  are  always  spherical  and  do  not  change  in  size. 

(3)  No  gravitational  force  acta  on  the  droplets. 

(4)  The  ratio  of  the  local  viscosity  to  the  free-stream  viscosity 
1^0^  appearing  in  equation  (4),  is  unity.  This  aBSun5>tion  is  valid, 

because  viscosity  is  not  appreciably  affected  by  the  order  of  magnitude 
of  the  changes  in  temperature-found  along  a trajectory.  “ 


METHOD  OP  CALCULATION 

The  nonlinear  differential  equations  of  motion  (equations  (l))  are 
difficult  to  solve  by  ordinary  means,  because  tbe  velocity  components 
of  the  air  and  the  term  containing  the  coefficient  of  drag  are  variable 
at  all  points  in  the  flow  field  and  cannot  be  related  in  closed  form. 
Simultaneous  solutions  for  the  two  equations  were  obtained  with  a mechan- 
ical analog  bsised  on  the  principle  of  a differential  analyzer  (refer- 
ence 3).  The  answers  were  obtained  in  the  form  of  plots  of  the  droplet 
trajectories  with  respect  to  the  cylinder,  as  shown  in  figure  2. 
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The  second-quadrant  section  of  the  cylinder  Is  outlined  In  fig- 
ure 2.  As  Is  shovm  In  the  flgwe,  the  ordinate  scale  was  expanded 
approximately  four  times  the  abscissa  scale  in  the  mechanical  analog  in 
order  to  gain  accuracy  in  determining  the  points  of  impingement  of  the 
droplets  on  the  cylinder  surface.  The  droplets  were  assumed  to  start  at 
an  infinite  distance  ahead  of  the  cylinder  at  an  initial  ordinate  of 
yQ  (fig.  l),  but  only  the  last  5 radii  of  distance  in  the  x-direction 

were  plotted  by  the  machine.  The  starting  conditions  at  x = -5  were 
obtained  by  use  of  equations  described  in  references  2 and  3.  The  ref- 
erence angle  0 in  figure  2 is  measured  clockwise  from  the  negative 
x-axis  of  the  cylinder. 

Trajectories,  such  as  those  shown  in  figure  2,  computed  for  several 
combinations  of  values  of  K and  Rcq  permit  the  evaluation  of  the 

amount  and  extent  of  imter-droplet  ln^singement  on  cylinders.  The  cylinder 
data  are  presented  herein  in  terms  of  the  dimensionless  parameters  K 
and  cp  used  in  reference  3.  The  pgirameter  K (equation  (2))  has  been 
termed  the  inertia  parameter  because  its  magnitude  directly  reflects  the 
external  force  reqMred  on  a droplet  to  cause  a deviation  from  the 
original  line  of  motion  of  the  droplet.  The  dimensionless  parameter  cp 
is  defined 


The  magnitude  of  P is  a measure  of  the  deviation  from  Stokes ' law  for 
the  forces  acting  on  the  water  droplets.  Stokes’  law  was  derived  for 
slow  translatory  motion  of  a small  sphere  in  an  incompressible  vlsco\is 
fluid  and  applies  precisely  in  the  limiting  value  of  cp  = 0. 

The  parameters  K and  CP  together  are  sufficient  to  define  the 
conditions  of  a particular  impinging  droplet.  Because  of  the  nature  of 
the  dimensionless  parameters,  an  infinite  ntnnber  of  combinations  of 
values  of  the  variables,  such  as  free-stream  velocity,  cylinder  size, 
droplet  size,  and  others  (equations  (2),  (3),  and  (6)),  would  apply  to 
the  particular  value  of  the  dimensionless  parameter  under  consideration. 
A system  of  equations  for  the  evaluation  of  dimensionless  parameters  in 
terms  of  variables  with  units  commonly  engjloyed  in  aeronautics  is  pre- 
sented in  appendix  C of  reference  3. 

The  compressible -air  velocity  components  in  the  flow  field  (equa- 
tions (5))  were  calculated  at  a Mach  number  of  0.4,  because  the  effect 
of  compressibility  for  completely  subsonic  flow  is  considered  greatest 
at  or  near  the  flight  critical  Mach  number.  The  value  obtained  for  the 
flight  critical  Mach  n\miber  depends  on  the’n\imiber  of  correction  terms 
applied  to  the  incompressible  flow  field  suid  appears  to  approach  0.4  as 
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a limit.  This  value  of  Mach  number  correbi<ond8  to  a speed  of  approxi- 
mately 290  miles  per  hour  at  an  altitude  of  10^000  feet.  Results  pre- 
sen-fed  herein  involve  the  maximum  effect  of  compressibility,  which 
restricts  the  free-stream  velocity  to  a value  compatible  with  a Mach 
number  of  0.4. 

Because  of  the  dlff ic\iltles  in  making  precise  corrections  to  the 
coefficient-of-drag  term  in  equations  (l)  and  to  the  density  terms  in 
equation  (4),  the  problem  was  solved  in  two  steps.  The  first  step  was 
to  survey  a range  of  values  of  the  parameters  K and  cp  to  determine 
the  conditions  most  affected  by  compressibility.  For  this  survey  the 
term  Cj^e/24,  containing  the  coefficient  of  drag  for  the  droplets, 

reqxiired  in  equations  (l),  was  assumed  to  be  a function  of  the  local 
Reynolds  number  Re  only  and  not  related  to  the  Mach  number.  The  valties 
for  Cj)Re/24  were  obtained  from  tables  in  reference  2.  Also,  the  den- 
sity ratio  appearing  in  equation  (4)  was  assumed  to  be  unity  for  this 
s\rrvey.  The  second  step  was  to  apply  first-order  corrections  to  the 
coefficient-of-drag  term  and  to  equation  (4)  for  density  variations 
along  the  particular  trajectory  that  showed  the  largest  effect  caused 
by  the  compressible-air  flow  field.  The  difference  in  the  trajectories 
caused  by  the  correction  was  less  than  the  probable  error  of  the  amalog. 


RESULTS  AMD  DISCUSSION 
Fi’esentatlon  of  Results 

The  solid  lines  in  figure  2 are  the  trajectories  for  droplets  in  a 
compressible -air  flow  field  calculated  at  the  flight  critical  Mach  number 
for  the  cylinder  (Mq  = 0.4),  and  the  dashed  lines  are  the  trajectories 

for  droplets  in  an  incompressible  flow  field.  The  set  of  lines  labeled 
A are  trajectories  of  droplets  with  large  momentum,  that  is,  large  in 
size  and  moving  with  high  speed.  The  momentum  of  these  droplets  is  so 
large  that  the  trajectories  are  nearly  straight  lines.  One  combination 
of  operating  conditions  for  \rtiich  the  trajectories  labeled  A would  be 
applicable  Involves  a cylinder  l/2  inch  in  diameter  moving  with  a free- 
stream  velocity  of  290  miles  per  hour  through  a cloud  composed  of  uniform 
droplets  approximately  26  miicrons  in  diameter  at  a 10,000-foot  edtitude. 
The  trajectories  labeled  A are  tangent  to  the  cylinder  at  the  point  of 
interception.  All  droplets  with  the  same  values  of  parameters  K and  p 
having  trajectories  below  these  trajectories  strike  the  cylinder; 
whereas,  ell  droplets  with  the  same  values  of  p and  K having  trajec- 
tories above  these  lines  will  miss  the  cylinder.  The  initial  ordinate 
y^  of  the  droplets  at  x = - » (fig.  1)  is  given  in  figure  2. 

A comparison  of  the  trajectories  computed  for  the  compressible-air 
flow  field  with  those  for  an  incompressible  flow  field  shows  that  for 


droplets  with  large  values  of  the  inertia  parameter  K (trajectories  A, 
fig.  2)  the  effect  of  compressibility  is  inperceptible.  The  effect  of 
compressibility  is  also  negligible  on  trajectories  of  droplets  with  ^ow 
values  of  K,  as  shown  by  the  set  of  tangent  trajectories  labeled  B. 

These  trajectories  would  apply,  for  example,  to  a cylinder  1 inch  in 
diameter  moving  with  a free-stream  velocity  of  290  miles  per  hoxir  through 
a cloud  composed  of  uniform  droplets  approximately  4 microns  in  diameter 
at  a 10,000-foot  altitude. 

The  greatest  effect  of  conpressibility  on  tangent  trajectories  is 
found  in  the  middle  range  of  the  droplet  inertia  parameter  K,  as  shown 
by  the  set  of  trajectories  labeled  C in  figure  2.  Both  trajectories 
are  tangent  to  the  cylinder,  but  the  droplets  had  to  start  at  different 
ordinates  in  order  to  arrive  tangent.  For  the  droplets  having  trajec- 
tories labeled  C,  the  collection  efficiency  (Ejj  = yo^m)>  defined  as  the 

ratio  of  the  actusJ.  water  impinging  on  the  cylinder  to  the  total  water 
in  the  vol\jme  swept  out  of  its  path  by  the  cylinder,  is  0.501  if  com- 
pressibility is  accounted  for,  as  compared  with  0.514  when  compressibility 
is  ignored. 


The  collection  efficiency  Em  for  a cylinder  in  a compressible-air 
flow  field  at  approximately  the  flight  critical  Mach  numiber  (Mq  = 0.4) 

is  presented  in  figure  3 as  a fimctlon  of  inertia  parameter  K and 
dimensionless  parameter  «p. 


The  effect  of  compressibility  on  the  collection  efficiency  is  shown 
in  figure  4 as  a percent  difference.  The  difference  is  obtained  from 
the  relation 


percent  difference 


^,i  " ®m. 


(100) 


j 

where  the  values  for  are  obtained  from  the  collection  efficiency  i 

presented  in  reference  3 for  the  incompressible  flow  field.  The  greatest 
difference  is  about  3 percent  at  the  conditions  of  K = 5 and 
Cp  = 50,000.  The  following  tabxilation  summarizes  the  effect  of  compres- 
sibility as  a percent  difference  between  results  computed  in  a compres-  i 

sible  and  an  incompressible  flow  field.  Typical  values  of  dimensions 
and  fli^t  conditions  are  Involved  in  the  examples  given.  All  calcula- 
tions were  made  for  290  miles  per  hour,  which  is  near  the  critical  speed 
for  cylinders  at  the  10,000-foot  RACA  standard  altitude. 
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Inertia 

parameter 

K 

Cylinder 

diameter 

(in.) 

Droplet 

diameter 

d 

(microns) 

Dimension- 

less 

parameter 

cp 

^,i  ” ^m,c  , . 

— u.  (100) 

%,c 

(percent) 

0.5 

1/2 

3 

740 

0.3 

.5 

1 

4 

1,480 

.4 

.5 

3 

7 

4,440 

.5 

.5 

6 

9 

8,880 

.5 

~1 

1/2 

4 

740 

.6 

1 

1 

5 

1,480 

.8 

1 

3 

9 

4,440 

.9 

1 

6 

13 

8,880 

1.0 

5 

1/2 

9 

740 

1.5 

5 

1 

12 

1,480 

1.9 

5 

3 

21 

4,440 

2.1 

5 

6 

30 

8,880 

2.3 

5 

18 

50 

27,650 

2.8 

5 

30 

70 

44,400 

2.9 

20 

1/2 

17 

740 

.5 

20 

1 

24 

1,480 

.6 

20 

3 

40 

4,440 

.7 

20  ' 

6 

60 

8,880 

.8 

Dlscxisslon  of  Results  j 

Prom  the  results  of  the  tabvilatlon  and  of  figure  4,  It  may  be  con-  j 

eluded  that  for  a cylinder  the  effect  of  coD5)ressibility  of  the  air  on  i 

the  droplet  trajectories  is  not  very  great  v^)  to  the  flight  critical 
Mach  number;  and  for  most  practical  problems  involving  measurements  of 
droplet  impingement  on  cylinders,  the  con^jressibility  effect  may  be 
ignored. 

The  reason  that  the  effect  of  conpressibility  on  droplet  trajec — j 

tories  is  nearly  negligible  is  foiind  in  the  conqjarison  of  the  j 

i compressible-air  flow  field  with  the  incompressible  flow  field.  The  j 

effect  of  compressibility  is  roost- pronovuiced  on  those  streamlines  close  ; 

i to  the  cylinder.  The  compressibility  effect  decreases  very  rapidly  with 

I distance  away  from  the  surface  of  the  cylinder  and  everywhere  becomes  ! 

negligible  1 cylinder  radius  away  from  the  surface.  Furthennore,  the 
effect  of  compressibility  is  of  consequence  only  over  the  top  portion  of 
the  cylinder,  as  shown  in  figure  5.  The  comparison  between  the  compres- 
sible and  incompressible  tangential  air -flow  velocities  on  the  cylinder  I 

surface  (fig.  5)  shows  that  the  effect  is  pronounced  only  for  6 >60°.  j 

The  angle  is  measured  from  the  same  reference  line  as  in  figure  1. 

i 

I • I 
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Droplets  forming  the  solid-line  trajectory  labeled  A in  figure  2 
are  acted  on,  during  the  last  portion  of  the  trajectory,  by  air  stream- 
lines that  have  a different  pattern  from  the  air  streamlines  for  an 
incompressible  flow  field;  however,  the  inertia  of  these  droplets  is  so 
large  that  the  last  portion  of  their  path  is  not  Influenced  appreciably 
by  the  con5»ressibility  of  the  air  flow.  The  droplets  with  low  inertia 
(trajectories  labeled  B,  fig.  2)  are  always  in  the  flow-field  region 
where  there  is  no  appreciable  difference  between  the  con5>ressible  and 
incompressible  air  flows.  The  droplets  in  the  middle  range  of  inertia 
have  the  last  portion  of  the  trajectories  in  a region  where  the  effect 
of  conpresslbility  on  the  air  flow  is  of  some  consequence.  The  Inertia 
of  these  droplets  is  low  enough  to  be  influenced.  Droplets  with  trajec- 
tories labeled  D and  E have  the  same  middle  range  of  Inertia  as 
those  with  trajectories  C.  Trajectories  D and  E are  not  Influenced 
by  compressibility  for  the  same  reason  that  trajectory  B was  not 
influenced. 

The  sample  trajectories  presented  in  figure  2 sianmarlze  the  effect 
of  comipresslbility  over  the  complete  range  of  values  studied  for  the 
parameters  K euid  *P.  The  small  effect  of  compressibility  is  most 
pronounced  on  collection  efficiency  (figs.  3 and  4).  The  effect  on  the 
maxlmnmn  angle  of  impingement  was  within  the  accuracy  of  the  analog. 

The  effect  on  trajectories  intermediate  between  the  tangent  trajectory 
and  the  x-axis  is  imperceptible  (trajectories  D and  E,  fig.  2).  For 
this  reason,  the  values  of  the  local  rate  of  impingement  (the  term  p * I 

in  reference  3)  are  not  affected  by  compressibility.  I 


Extension  of  Cylinder  Results  to  Airfoils 

The  extension  of  the  results  obtained  with  a cylinder  to  an  airfoil 
is  Justifiable,  because  the  incompressible  flow  fields  of  cylinders  and 
airfoils  are  similarly  altered  by  compressibility,  A comparison  between 
the  compressible  and  incompressible  velocities  on  the  surfaces  of  two 
symmetrical  Joukowski  airfoils  is  shown  in  figure  6.  The  ctnrves  of  fig- 
ure 6 were  obtained  from  results  presented  in  reference  8.  The  ratio 
of  the  surface  velocity  to  the  free -stream  velocity  is  plotted  as  a 
function  of  the  chordwise  station,  given  as  a ratio  to  the  chord  length. 
The  results  of  figure  6 indicate  that  the  effect  of  compressibility 
becomes  less  pronounced  as  the  thickness  ratio  t of  the  airfoil  is 
decreased  and  the  critical  Mach  number  is  increased. 

As  was  found  for  the  cylinder,  the  greatest  effect  of  compressibility 
is  in  the  vicinity  of  the  thickest  part  of  the  airfoil.  Analysis  based 
on  the  work  presented  in  references  5 and  6 indicates  that  under  normal 
icing  conditions  the  impingement  occurs  on  the  nose  section  of  the  air- 
foil. As  with  the  cylinder,  the  nose  section  is  not  affected  to  a sig- 
nificant degree  by  compressibility.  Those  droplets  that  impinge  on  the 
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thick  part  of  the  airfoil  have  such  large  inertia  that  small  differences 
in  air-flow  velocity  components  close  to  the  airfoil  do  not  appreciably 
affect  their  trajectories. 

Aerodynamic  considerations  show  that  the  superposition  of  circulation 
on  the  air-flow  field,  as  would  be  caused  by  placing  the  airfoil  at  an 
angle  of  attack,  would  tend  to  cause  a shift  in  the  location  of  the  max- 
imum difference  between  the  compressible  and  incompressible  air  flows. 

The  circulation  would  shift  the  maximum  difference  toward  the  nose  on 
the  top  sxurface  and  toward  the  tail  on  the  bottom  surface,  because  cir- 
culation-increases the  air  velocity  around  the  top  surface,  especially 
near  the  nose,  auad  decreases  the  velocity  near  the  bottom  surface. 

However,  a change  in  angle  of  attack  also  causes  a shift  in  the  area  of 
impingement  in  the  same  direction,  that  is,  toward  the  nose  on  the  top 
surface  and  toward  the  tall  on  the  bottom  surface;  therefore,  the  area  of 
impingement  would  probably  remain  in  the  region  where  the  effect  of  com- 
pressibility is  of  no  practical  consequence. 

The  curves  of  figxire  6 are  for  (Joukowskl  airfoils,  which  are  more 
blunt-nosed  than  modern  high-speed  airfoils.  The  chordwise  location  of 
the  maximum  velocity  occurs  at  approximately  x = 0.1  for  the  Joukowski 
airfoil;  whereas,  the  chordwise  location  of  the  maximum  velocity  occurs 
from  x=0.3  to  x=0.5  on  high-speed  airfoils,  depending  on  the 
design  lift  coefficient.  The  area  of  droplet  iiip>ingement , however,  is 
at  the  nose  section  for  a high-speed  airfoil  (reference  6),  as  in  the 
case  of  the  Joiokowskl  airfoil  (reference  5).  The  arguments  used  in  con- 
clxidlng  that  the  effect  of  the  conpressibillty  of  air  is  negligible  on 
droplet  trajectories  inpinging  on  a Jovikowski  airfoil,  for  which  the 
flow  field  is  a direct  transformation  from  a cylinder,  are  applicable  to 
modern  high-speed  airfoils,  because  the  air-Tlow  regions  with  the  great- 
est conpressibility  of  the  air  occur  on  the  airfoil  well  back  of  the 
regions  of  droplet  impingement. 


CONCLUSIONS 

The  following  conclusions  are  based  on  the  calculation  of  droplet 
trajectories  around  a cylinder  in  both  an  inconpressible  flow  field  Euid 
a compressible -air  flow  field  at  the  flight  critical  Mach  number  and  on 
a study  of  the  air  velocities  around  airfoils; 

1.  For  a cylinder  the  effect  of  compressibility  of  the  air  on  the 
droplet  trajectories  is  not  very  great  up  to  the  flight  critical  Mach 
number,  smd  for  most  practical  problems  involving  measurements  of 
droplet  impingement  on  cylinders  the  conpressibility  effect  may  be 
ignored. 
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2.  The  fact  that  the  effect  of  con5)res8lhility  on  the  impinging 
droplet  trajectories  is  negligible  also  applies  to  airfoils. 


lewis  Plight  Propiilslon  Laboratory 

National  Advisory  Committee  for  Aeronautics 
Cleveland,  Ohio,  September  11,  1952 
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Cylinder  angle,  6 


Flgttre  5.  - Contpeo'lsoD  of  surface  velocity  on  a cylinder 
between  incompressible  flow  field  and  compressible-air 
flow  field  at  approximately  flight  critical  Mach 
numiber. 
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Tangential  surface  velocity,  Ug,  ratio  to  free-stream  velocity 
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Figure  6.  - Cong>arlBon  of  surface  velocities  on  symetrical  Joukovski 
airfoils.  Angle  of  attack,  O*’. 


HACA-l«ii(l«r  • 3.6-M  - 1000 


f,; 


